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1.  INTRODUCTION 


Hie  dectrothennal-chemical  (ETQ  gun,  generically  shown  in  Figure  1 ,  is  a  propulsion  concept  which 
utilizes  a  low-mass,  high-energy  plasma  to  initiate  and,  hopefully,  control  the  combustioiWaponzation  of 
the  woiking  fluid  (propellam)  during  the  ballistic  cycle.  Controlling  combustion  (exothermic  working 
fluid)  or  vaporization  (endothermic  woiking  fluid)  of  the  woiking  fluid  in  the  ETC  gun  is  necessary  in 
order  to  tailor  die  pressure-time  profile  in  the  guiL  Theoretically,  tailoring  die  pressure-time  profile  to 
obtain  a  "flatter”  and  extended  pressure  curve  should  result  in  enhanced  performance  (increased  muzzle 
energy)  and  provide  a  "softer"  launch  environmem  to  enhance  projectile  integrity,  especially  for  "smait" 
projectiles.  In  addition,  precise  control  of  the  pressure  history  should  also  allow  for  the  required  projectile 
velocity  repeatability  required  of  indirect  fire  support  (artillery)  applications.  An  essential  facet  of 
understanding  the  control  of  the  process  throu^  the  plasma-propellant  interaction  is  an  accurate 
description  of  the  energy  release  rate  (J/s). 

In  previous  worit  (Wren  and  Obeile  1990,  1992),  the  authors  derived  and  implemented  an  inverse 
analysis  of  experimental  ETC  gun  firing  data  to  determine  if  a  relationship  between  electrical  energy  input 
and  energy  release  rates  was  indicated.  Fbr  that  analysis,  variable  thermochemical  properties  of  the 
propellant  gas/j>lasma  mixture  were  assumed  to  be  dependent  only  on  the  constitutive  components  of  the 
propellam  and  the  electrical  energy  density  (kJ/g;  ratio  of  total  electrical  energy  input  to  consumed 
propellant).  Implementation  of  this  variable  theimodiemistiy  in  the  computer  model  was  through  the  use 
of  "look-up"  uUes,  similar  to  Table  1  for  JA2  propellant  (Bunte  and  Oberie  1989).  Values  in  the  tables 
were  obtained  using  the  thennodynamic  equilibrium  code  BLAKE  (Freedman  1982)  at  a  fixed  loading 
density  (Id)  (g/cm^;  ratio  of  consumed  propellant  to  flee  volume)  of  0.2  g/cm^. 

However,  a  detailed  analysis  of  the  results  from  the  study  identified  several  potential  shortcomings 
related  to  the  use  of  look-up  tables  for  the  theimochemical  properties.  First,  die  values  of  the 
theimochemical  properties  are  dependent  on  the  loading  density,  as  illustrated  in  Table  2  for  JA2 
propellant 

Fortunately,  it  has  been  shown  (Robbins  1991)  that  for  traditional  solid  propellam  modeling, 
consideration  of  the  loading  density  dependence  of  values  of  the  thermochemical  properties  has  no 
significant  (<1%)  impact  on  simulation  results,  even  through  loading  density  (ratio  of  mass  of  gas  to  gas 
volume  in  g/cm^)  can  vary  over  a  large  range,  as  illustrated  in  Figure  2.  Similar  results  have  yet  to  be 
demonstrated  for  ETC  simulations.  In  fact,  for  the  ETC  gun,  the  problem  of  determining  the  energy 
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Figure  1.  Schematic  of  ETC  gun. 


Table  1.  Variable  Thennochemical  ftopcrties  for  JA2  as  a  Function  of 
Electrical  Energy  Density 


Electrical 
Energy  Input 
(kJ/g) 

Temperature 

(K) 

Impetus 

(J/g) 

Molecular 

Weight 

3,424 

1,144 

24.886 

1 

3,959 

1,334 

24.686 

2 

4,401 

1,504 

24.337 

.  3 

4,776 

1,663 

23.882 

4 

5.113 

1,819 

23.378 

5 

5,429 

1,975 

22.857 

Specific 

Energy 

(J/g) 


1.2254 

1.2219 

1.2219 


5,075 

6,012 

6,778 

7.421 

7,996 

8,528 


1.2241 

1.2275 

1.2316 


Tabled  Loading  Doisity  Dqjendence  of  Values  for  Tbennocfaeinical  Properties 
Computed  by  BLAKE  for  IA2  Propellam 


T  /Mrfing 


.1000 


0.3000 

0.4000 


0.6000 


Teiiq> 

(K) 


3371 

3395 

3,406 

3.417 

3.425 

3.431 


Inqretus 

(J/g) 


1.148.3 


Molecular 

Weight 

Gas 


24.704 

24.749 

24.664 

24.797 

24.820 

24.846 


Co-Vol 

(an^/g) 

Frozen 

Y 

Spedfic 

Energy 

(f/g) 

1.031 

13229 

5.090 

0.991 

13254 

5,059 

0.948 

13305 

4.962 

0.903 

13379 

4.816 

0.858 

13474 

4.367 

0.814 

13589 

4.435 

JA2,  120mm 
Constont  Thermochems 
at  looding  density 
of  0.2  g/cc 
Averoge  Id  =  0.236 


Time  (ms) 


JJii'iiMilulLLtil’l'ii 


II  I 


release  rate  is  further  complicated  by  the  staged  introduction  of  electrical  energy  which  produces 
combustion  gases  with  even  larger  changes  in  thermochemical  values  compared  to  solid  propellants  alone 
(see  Tables  1  and  2).  Thus,  it  appears  that  for  ETC  interior  ballistic  GB)  calculations,  both  the  loading 
and  electrical  energy  densities  should  be  considered  in  determining  the  energy  release  rate.  Although 
some  earlier  (Oberle  1989;  Gough  1989;  Wren  and  Oberle  1990. 1992)  investigations  of  ETC  performance 
considered  variations  in  thermochemical  values  (energy  release  rate)  as  a  function  of  electrical  energy 
density,  no  work  has  incorporated  both  variations  in  loading  density  and  electrical  energy  density. 

A  second  and  potentially  more  serious  shortcoming  to  the  use  of  a  lo(dc-up  table  is  the  accuracy  of 
the  table  values  as  the  electrical  energy  density  increases  to  produce  average  gas  temperatures  above 
5.(XX)  K.  The  thermochemical  values  typically  utilized  in  IB  codes  are  flame  temperature,  chemical 
energy,  and  covolume  based  on  polynomial  fits  to  experimental  values  of  specific  heat  at  constant 
pressure,  C^.  In  the  case  of  BLAKE,  the  JANAF  data  for  for  all  of  the  produa  gases  are  computed 
using  statistical  mechanics  based  on  spectroscopic  data.  Although  these  calculations  can  be  readily 
extended  well  above  5,000  K,  tables  of  these  data  traditionally  terminate  aroimd  5,000  K.  Thus,  the 
fittings  of  the  data  are  valid  only  to  5.(X)0  K.  The  polynomial  fitting  fimctions  for  take  on  values  to 
plus  or  minus  infinity  outside  the  fitting  range;  will  behave  similarly.  Tlieir  ratio  y  will  approach  one 
in  the  limit,  and  ttie  value  of  y  above  5,000  K  is  thus  useless  for  the  ETC  a^lication  (Freedman  1991), 
since  gas  temperatures  for  ETC  (xopellants  may  be  high,  at  least  locally,  where  the  plasma  temperature 
is  typically  in  the  range  of  10,(X)0-15,(X)0  K.  Thermochemical  values  for  high  electrical  energy  densities, 
and  hence,  high  temperatures,  are  thus  ill-defined  and  based  on  extrapolated  polynomial  fits. 

Traditionally,  the  chemical  energy  is  determined  primarily  from  impetus  and  the  ratio  of  specific  heats 
(Y).  Since  y  is  directly  related  to  specific  heat  at  constant  pressure  and  is  supplied  to  BLAKE  as 
experimental  data,  it  might  be  expected  that  extrapolated  values  of  y  are  particularly  poorly  behaved.  In 
the  previous  effort  (Wren  and  Oberle  1990,  1992)  to  determine  energy  release  rates,  an  adequate 
description  of  y  above  5,000  K  was  a  major  shortfall.  In  addition,  the  y  calculated  by  BLAKE  is  for  the 
gas  phase  only  and  does  not  address  the  use  of  solid  particles  used  in  some  experimental  ETC  propellants. 
Previous  IB  calculations  have  shown  that  traditional  measures  of  "goodness"  of  the  propellant  may  not 
be  ^licable  to  ETC  propellants  containing  large  amounts  of  solid  particles  such  as  aluminum  (White  and 
Oberle  1989). 
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Therefore,  to  address  the  question  as  to  whether  fixed  thermochemistiy  will  be  sufficient  to  describe 
the  ETC  IB  process  and  bypass  the  use  of  7.  the  thetmochemical  code  BLAKE  has  been  directly  linked 
to  a  lumped  parameter  IB  model  and  an  inverse  model  developed  eariier.  The  computed  value  of  the 
space-mean  pressure  is  sui^lied  directly  to  the  IB  and  inverse  codes,  thereby  directly  including  the  effect 
of  the  loading  density  and  electrical  energy  density  and  obviating  the  need  for  impetus  and  y  in  the 
caictilarinn  of  energy.  Although  the  space-mean  pressure  is  also  related  to  the  experimental  values  of 
specific  heat  at  constant  pressure,  it  is  expected  to  be  much  more  reliable  since  it  is  a  Amction  of  several 
well-defined  variables  (Freedman  1991). 

2.  DESCRIPTION  OF  THE  MODELS 

The  thetmochemical  code  BLAKE  was  revised  to  serve  as  a  subroutine  to  a  main  calling  program. 
BLAKE  requires  the  constituem  data  as  well  as  values  of  electrical  energy  density  (electrical  energy  irgait 
in  kilojoulesAnass  of  propellant  consumed  in  grams)  and  loading  density.  The  thermochemical  code  then 
returns  the  space-mean  pressure  of  the  gas  assuming  equilibrium  thermodynamics. 

Two  implementations  of  the  thetmochemical  code  are  explored  in  this  rqoort  The  first  is  a  traditional 
time-marching  IB  formulation  using  a  Lagrange  gradient  and  expressed  in  a  form  to  utilize  the  space-mean 
pressure,  denoted  as  IBBLAKE.  This  formulation  uses  as  irgwt  the  gun  geometry,  projectile  mass  and 
resistive  pressure  profile,  propellant  geometry,  propellant  tiierroochemistry,  and  electrical  energy  history. 
By  iterating  on  the  mass  of  propellant  consumed  in  a  time  step,  which  is  passed  to  BLAKE  along  with 
the  electrical  energy,  the  conservation  equations  are  satisfied.  Since  the  equilibrium  state  of  the  gas  at  any 
time  step  must  consider  the  woric  during  expansion  in  terms  of  projectile  kinetic  energy,  fluid  kinetic 
energy,  and  other  losses  (assumed  0.0  for  the  purposes  of  diis  study),  the  electrical  energy  is  reduced  by 
the  predicted  kinetic  energies  of  the  gas  and  projectile  on  any  time  step.  (Note:  This  may  result  in  a 
negative  electrical  energy  input  to  the  BLAKE  calculation.)  The  mass  of  prc^llant  consumed  is  then 
determined  and  the  pressures,  projectile  motion,  and  other  IB  information  are  ouqxit 

Secondly,  the  BLAKE  subroutine  was  integrated  into  an  inverse  code  (Wren  and  Oberle  1990, 1992) 
which  utilizes  experimental  data  in  order  to  infer  the  mass  of  propellant  consumed  at  any  time  step.  The 
linked  model  is  denoted  as  INVBLAK.  The  inverse  model  uses  as  input  experimental  values  of  projectile 
position,  electrical  energy  input  to  the  gun.  and  chamber  pressure.  Assuming  a  Lagrange  gradient,  it  is 
then  p(»sible  to  infer  ^ace-mean  pressure  in  the  experiment  since  the  gas  volume,  projectile  kinetic 


energy,  and  fluid  kinetic  energy  are  known.  The  BLAKE  subroutine  is  used  to  supply  a  space-mean 
pressure  based  on  an  electrical  energy  density  and  loading  density  for  the  propellant  gases.  By  iterating 
on  the  mass  of  propellant  consumed,  the  space-mean  pressure  in  the  experiment  and  the  space-mean 
pressure  predicted  by  BLAKE  are  simultaneously  satisfied.  The  ou^ut  is  then  the  mass  of  propellant 
consumed  at  time  step. 

3.  COMPARISON  WITH  SOLID  PROPELLANT  MODEL 

In  order  to  assess  the  effect  of  incorporating  updated  thermochemical  values  in  a  solid  propellant  only 
simiiiatinn,  die  gun  described  in  Table  3  was  modeled  with  IBHVG2  and  IBBLAKE.  To  perform  the 
simulation,  a  choice  concerning  the  bum  rate  and,  hence,  the  mass  generation  rate  for  the  solid  propellant, 
must  be  made.  For  both  IB  codes,  the  bum  rate  is  provided  in  the  form  r  =  hP".  If  fixed  thermochemical 
values  are  utilized  in  die  IB  code,  then  the  choice  for  the  bum  rate  is  straightforward.  Use  the  coefficient 
b  and  expcment  n  determined  via  closed  chamber  data  and  the  fixed  diermochemical  values.  (Note  that 
the  bum  rate  is  thus  a  function  of  the  thermochemical  values).  However,  if  variable  thermochemical 
values  are  to  be  used  in  the  IB  code,  then  the  bum  rate  should  be  adjusted  to  be  consistent  with  the  closed 
chamber  data  and  thermochemical  values.  Since  for  this  report  the  thermochemical  values  are  a  function 
of  loading  density,  it  was  necessary  to  adjust  the  tnim  rate  to  account  for  the  variability  in  Id  shown  in 
Figure  2.  One  approach  would  be  to  perform  an  extoisive  series  of  closed  chamber  firings  at  different 
loading  densities.  However,  such  a  series  of  firings  was  not  feasible  for  this  report 


Table  3.  Gun  Parameters  Used  in  Simulation 


Bore  Diameter.  120  mm 

Projectile  Mass:  11.4  kg 

Chamber  Volume:  9,750  cm^ 

Charge  Mass:  8.8  kg 

Projectile  Travel:  475  cm 

Propellant:  JA2 

Thus,  the  bum  rate  as  a  function  of  loading  density  is  estimated  by  the  following  approach.  The  final 
equation  for  computing  the  bum  rate  r  from  closed  chamber  data  is 


ih 

r  « - - 

pA 


(1) 
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where  ih  is  the  time  rate  of  change  of  mass,  p  the  material  density,  and  A  the  reacting  surface  area. 
Assuming  constant  propellant  density,  p,  and  a  neutral  (constant)  reacting  surface  area.  A,  the  bum  rate 
is  directly  proportional  to  the  time  rate  of  change  of  the  mass,  in. 


r  «  km.  (2) 

Now  the  ina»y  history  (mass  vs.  time)  is  determined  from  the  closed  chamber  pressure  history  by 
dfttftnnining  the  total  Chemical  energy  necessary  to  produce  the  observed  pressure  and  the  relation. 


total  chemical  energy  (TCE) »  mass  propellam  consumed  *  propeUam  specific  energy 

=  m  *  e  .  (3) 

If  variable  thermochemical  values  are  utilized,  then  the  propellant  specific  energy,  e.  wiU  also  be  a 
function  of  loading  density.  However,  the  total  diemical  energy  required  is  depeivlent  only  on  the 
observed  dosed  chamber  pressure. 

Thus, 

/?Iyy(/*j)  •  ey^(/* |)  «  TCE(P |) 

•  my2(/*i)  •  eu2(Pi)  ,  (4) 


vdiere  represents  an  observed  pressure,  and  ld\  and  Idl  values  are  associated  with  two  different  loading 
densities.  To  compute  M,  consider  two  pressures,  and  ^2*  measured  in  a  time  step  Ar.  From 
Equation  4, 


At 


At 


(5) 
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However,  the  specific  energy  is  a  function  of  loading  density,  not  pressure;  thus, 


or 


^ld\  ^ld\  •  ^Id2^ld2 


(6) 


'^Idl  ^id2 


Ccmibining  Equations  2  and  7, 


'ld\  ^^ld\  ^ld2 

•  _  » 

rid2  ^Idl 


(8) 


or 


_ ^ _ ^/£i _ 

r{ld  «  Q2glcm^)  e(ld  ■  0,2g/cm^) 


(9) 


assuming  that  the  closed  chamber  data  has  been  fitted  with  a  ime-part  bum  rate,  which  is  usually  the  case. 
Equation  9  is  used  to  adjust  the  bum  rate  as  a  function  of  loading  density.  For  the  single-perforated  grain 
under  consideration,  die  neutral  surface  area  assumption  is  felt  to  be  reasoiudile. 

The  results  are  shown  in  Table  4.  It  can  be  seen  diat  the  results  are  in  good  agreement,  indicating 
that  the  approximation  of  the  thermochemistry  of  the  propelling  gas  using  constant  thermochemical  values 
at  a  Id  of  0.2  g/crn^  and  a  pressure-dqiendent  bum  rate  law  based  on  a  Id  of  0.2  g/cm^  used  in  IBHVG2 
are  reasonable.  Thus,  it  appears  that  constam  thermochemical  values  at  a  Id  of  0.2  g/cm^  provides  a  good 
simulation  of  solid  prcqiellant  performance. 


TaUe  4.  Ccunparison  of  IBHVG2  and  IBBLAKE  for  Simulation  of  120-mm 
Gun  Using  Single-Perforated  JA2  Prt^Uant 


Maximum 

Pressure 

Muzzle  Velocity 

(MPa) 

(m/s) 

IBHVG2 

1,403 

IBBLAKE 

1,415 
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In  the  case  of  ETC  propellants,  energy  is  not  only  a  function  of  loading  density  but  a  function  of 
electrical  energy  density  as  well.  Figure  3  shows  a  gr^h  of  energy  vs.  Id  vs.  electrical  energy  density 
for  a  representative  candidate  ETC  propellant  of  80%  by  mass  of  HAN  and  methanol  in  a  stoichiometric 
mixture  and  20%  water.  Current  lumped  parameter  ETC  IB  models  using  variable  thermochemistry  utilize 
tables  of  energy  vs.  electrical  energy  density  at  a  constant  Id  of  0.2  g/cm^.  Figure  3  shows  significant 
variation  of  energy  of  up  to  30%  at  a  constant  td  -  0.2  g/cm^,  with  electrical  energy  densities  in  the  range 
of  0.4  kJ/g  to  2.0  kJ/g.  Thus,  it  might  be  expected  that  greater  differences  between  IBHVG2  and 
IBBLAKE  would  be  observed. 

In  order  to  determine  a  qualitative  difference  between  the  two  types  of  simulations  (and  since  a 
burning  rate  is  not  known  for  the  ETC  propellant  considered),  the  simulation  is  performed  with  all  input 
parameters  identical  to  those  used  previously  for  the  solid  propellam  except  that  the  theimochemical 
properties  are  changed  to  those  of  the  HAN,  methanol,  and  water  mixnire  and  the  burning  rate  is  not 
adjusted.  A  comparison  of  the  maximum  pressure  and  muzzle  velocity  is  shown  in  Table  S.  The 
maximum  breech  pressures  are  comparable,  and  the  projectile  velocity  differs  by  4%,  a  larger  difference 
than  noted  in  the  JA2  simulatiorL 

Table  S.  Comparison  of  IBHVG2  and  IBBLAKE  for  Simulation  of  120-mm 
Gun  Using  a  HAN/Methanol/Water  Mixture 


Maximum 

Pressure 

Muzzle  Velocity 

1 

(MPa) 

(m/s) 

1  IBHVG2 

350 

1,128 

1^  IBBLAKE 

349 

1,172 

Thus,  it  tqrpears  that  IB  models  using  constant  thermochemical  values  at  a  loading  density  of 
0.2  g/cm^,  combined  with  pressure-dependent  experimental  data  reduced  at  the  same  Id,  provide  good 
simulation  of  traditional  solid  propellant  guns  based  on  this  limited  study.  However,  ETC  propellants 
combined  with  electrical  energy  may  display  thermochemical  properties  which  are  more  strongly  affected 
by  loading  density  considerations.  Thus,  it  appears  that  traditional  IB  methods  of  approximating  the 
properties  of  the  propelling  gas  may  not  be  adequate  in  ETC  simulations. 
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^=^ms\ 
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if  HAN  and  Methanol  in  a  Stoichiometric 


4.  COMPARISON  WITH  INVERSE  MODEL 

An  inverse  model  was  developed  previously  (Wren  and  Obeile  1990,  1992)  using  variable 
thermochemisoy  at  Id  =  0.2  g/cm^  and  applied  to  a  number  of  experimental  shots.  The  model  requires 
ejqperimental  breech  pressure,  projectile  position,  and  electrical  energy  iigxit  at  a  time  step  and  infers  the 
mass  of  propellam  consumed  based  on  an  energy  balance  equation.  At  gun  pressures,  compression  of  the 
liquid  propellant  can  result  in  volume  changes  whidi  affect  the  calculation  of  consumed  mass.  Thus,  the 
model  has  been  extended  to  treat  compression  of  the  working  fluid  (propellant)  and  consider  losses  as  well 
as  to  extrapolate  dK  thermochemical  data  (Wren  arvl  Oberle  1992).  The  model  was  revised  for  this  study 
to  permit  the  use  of  BLAKE  as  a  subroutine.  The  mass  of  propellam  consumed  is  that  required  to  produce 
the  space>mean  (nessure  inferred  from  the  experimem. 

Electrical  energy  input  for  a  30-mm  experimetual  firing.  Shot  39,  performed  by  GT-Devices  (Gteig 
1990)  as  part  of  a  repeatal^ty  series  using  titanium  hydride  and  water  is  shown  in  Figure  4  by  the  dotted 
line.  The  {wevious  inverse  model  results  (Wren  and  Oberle  1992)  using  a  taUe  of  variable  thermochemical 
prt^rties  (see  TaUe  6)  at  a  fixed  Id  of  02  g/cm^  is  shown  in  Figure  4  by  the  solid  line.  In  previous  work, 
thermochemical  values  outside  die  range  of  TaUe  6  were  taken  to  be  the  first  and  last  values  in  the  table 
(Wren  and  Oberle  1990)  or  were  extrqrolated  (Wren  and  Oberie  1992).  However,  eariy  (during  the  first 
milliseoond)  in  die  IB  cyde,  electrical  energy  densities  are  predicted  to  be  quite  hi^  (above  10  kJ/g)  and 
temperatures  are  high  (above  S,(XX)  K).  However,  the  projectile  motion  is  also  not  well  described  during 
the  first  millisecond  in  Shot  39.  Hence,  the  estimate  of  mass  consumed  during  the  time  from  0.0  to 
1.0  ms  is  not  reliable.  It  is  noted  that  die  inverse  analysis  is  an  energy  balance  at  each  time  step 
independem  of  any  other  time  step.  Hence,  the  difficulty  with  prediction  of  mass  consumed  from  0.0  to 
1.0  ms  does  not  influence  die  piedicticm  at  a  later  time.  The  total  mass  of  propellam  in  the  experimem 
is  168.84  g.  As  indicated  by  die  solid  line  in  Figure  4,  approximatdy  150  g  of  propellam  is  inferred  to 
be  consumed  using  a  variaUe  thermochemistry  table  CTable  6). 

The  linked  inverse-BLAKE  code  (INVBLAK)  results  are  shown  by  the  triangles  in  Figure  4.  It  is 
seen  diat  significantly  less  fluid  is  inferred  to  be  consumed.  These  results  are  consistem  with  experimental 
observations  that  a  quantity  of  unbutned  fluid  was  present  at  the  conclusion  of  these  shots  (Greig  1990). 
For  comparison,  the  Id  vs.  time  required  by  BLAKE  to  match  the  experimental  space-mean  pressure  is 
shown  in  Figure  S.  It  is  noted  that  the  Id  stays  bdow  0.1  g/cm^,  a  regime  in  which  the  specific  energy 
is  higher  than  for  a  Id  of  0.2  g/cm^  (see  Figure  3).  The  experiment  has  a  large  unount  of  ullage  initially. 
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and  if  dl  tbe  propdlam  reacted,  the  final  Id  would  be  0.07  g/cm^  Thus,  the  inclusion  of  thennodiemical 
dependence  on  loading  and  electrical  eneigy  densities  in  die  inverse  modd  ^ipears  to  provide  a  better 
estimate  of  die  mass  of  propdlam  consumed  in  the  ETC  gun  firing  examined. 
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Table  6.  Variable  Thennodiemical  Values  for  50%  Ti3H2  and  50%  H2O  Plus  Electrical  Eneigy 
at  U  «  02  g/cm^ 


Electrical  Energy 
Input 
(U/g) 

Temperature 

(K) 

Impetus 

m 

Molecular 

Wdght 

Y 

Spedfic 

Eneigy 

(J/g) 

2 

2,162 

6792 

6.716 

12722 

2200 

3 

2.630 

826.7 

6.714 

12722 

3240 

4 

3,078 

969.0 

6.704 

12553 

3.970 

5 

3,502 

1.106.3 

6.682 

1.2443 

4.670 

6 

3,899 

12392 

6.645 

12370 

5230 

7 

4265 

1268.6 

6296 

12325 

5250 

8 

1.495.8 

6246 

12302 

6220 

9 

4,917 

1,622.4 

6.506 

12294 

10 

5210 

1.749.4 

6.487 

12313 

7^® 

Shot  39 


Noik  ExperiaeM  elMirical  energy  inint  (4k>i^  Man  of  {vopeOani  connimed  beMd  on  mvoae  model  with  uUe  of 
dieniiochemical  vahiei  (nlid).  Man  of  ptopettau  cgnsaaed  bend  on  invene  model  to  BLAKE  (tiinglet). 


Figure  4.  Inverse  modd  results  for  30-mm  titanium  hydride  and  water  firing. 
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Rgiife5.  Loxteu; <teisitY  w-  time  for  go-mm  tiaroiuiiL hydride  and  w«er firin£  based  on 
inverse  model  linked  to  BLAKE 


5.  CONCLUSIONS 


A  thennodKoiical  equilitniuin  code,  BLAKE,  has  been  directly  linked  to  both  an  IB  code  and  an 

♦  _ 

inveise  code  in  Older  to  assess  the  potential  inqxovement  in  describing  ttie  energy  rekare  nie  in  ETC  gun 

moddiqg.  The  inverse  analysis  attempts  to  detennine  the  decomporition  or  gas  generation  me  necessary 

to  satisfy  an  energy  balance  equation  based  on  experimental  data. 

Hie  analysis  was  apfdied  to  die  following:  1)  a  simulation  of  a  standard  120-imn  gun  widi  JA2 
propellant;  2)  a  simulation  of  a  standard  120-mm  gun  widi  a  potential  representative  ETC  propellant;  and 
3)  an  inverse  model  of  a  30-mm  eiqreiimental  firing  using  titanium  hydride  and  water  as  die  propellant 
(working  fluid).  The  results  suggest  die  following: 

(1)  The  diermodiemical  code  BLAKE  can  function  as  a  subroutine  to  an  IB  code. 
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(2)  Constant  values  of  thennochemical  properties  and  bum  rates  used  in  solid  propellant  modeling 
yield  comparable  results  to  the  linked  thennochemical-IB  model. 

(3)  There  is  a  greater  difference  between  an  IB  simulation  linked  to  BLAKE  of  an  ETC  propellant 
and  an  IB  simulation  using  constam  thermochemical  values  than  with  traditional  solid  propellants. 
The  muzzle  velocities  differ  by  4%.  a  larger  difference  than  the  solid  propellant  considered. 

(4)  A  consideration  of  variable  thermochemistry  via  a  direa  link  to  BLAKE  in  an  inverse  model  (as 
opposed  to  tables  of  theimochemical  values  based  on  a  constant  loading  density  of  02  g/cm^) 
makes  a  significant  difference  in  the  inferred  values  of  mass  of  propellant  consumed. 

(5)  Treatment  of  the  dependence  of  thermochemistry  on  both  electrical  energy  density  artd  loading 
density  appears  to  provide  a  better  model  of  the  state  of  the  propelling  gas  in  ETC  ai^ications. 
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Bldg.  3  IS 

Fort  Belvoir.  VA  22060-5606 

1  Commander 

U.S.  Army  TRAC  -  Fort  Lee 
Defense  Logistics  Studies 
F0rtLee.VA  23801-6140 

1  President 

U.S.  Army  Artillery  Board 
F0rtSai.OK  73503 

1  Commandant 

U.S.  Army  Command  and  General  Staff  College 
FOn  Leavenworth.  KS  66027-5200 

1  Commandant 

U.S.  Army  Special  Warfare  School 
ATTN:  Rev  and  Tng  Ut  Div 
FOrt  Bragg.  NC  28307 


No.  of 
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1  Commander 

Radford  Army  Ammunition  Plant 
ATTN:  SMCRA-QA/HI  Library 
Radford.  VA  24141 

2  Commander 

U.S.  Army  Foreign  Science  and  Technology  Center 
ATTN:  AMXST-MC-3. 

S.  LeBeau 

C.  Beiter 
220  Seventh  Sl,  NE 
Charlottesville,  VA  22901 

1  Commandant 

U.S.  Army  Field  Artillery  Center  and  School 
ATTN:  ATSF-(X)-MW,  B.  Willis 
FortSULOK  73503 

3  Deputy  Commander 
Strategic  Defense  Command 
ATm:  SFAE-SD-HVL, 

S.  Smith 
LTCKec 

D.  Uanos 
P.O.  Box  1500 
Huntsville.  AL  35887-3801 

1  Naval  Sea  Systems  Command 
Department  of  the  Navy 
ATTN:  CSEA,  CDR  Dampier 
06KR12 

Washington.  DC  20362-5101 

1  Office  of  Naval  Research 
ATTN:  Code  473,  R.  S.  Miller 
800  N.  Quincy  St. 

Arlington,  VA  22217 

2  Commander 

Naval  Sea  Systems  Command 
ATTN:  SEA62R 
SEA  64 

Washington.  DC  20362-5101 

1  Commander 

Naval  Air  Systems  Command 
ATTN:  AIR-954,  Technical  Library 
Washington,  DC  20360 
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Naval  Research  Laboratory 
Technical  Lilxary 
Washington.  DC  2037S 

Commander 

Naval  Surface  Warfare  Center 
ATTN:  J.  P.  Consaga 
C.  Gotzmer 

Silver  Spring,  MD  20902-5000 
Commander 

Naval  Surface  Warfare  Center 
ATTN:  K.  Kim.  Code  R-13 

R.  Bemecker.  Code  R-13 
Silver  Spring.  MD  20902-5000 

Conunander 

Naval  Surface  Warfare  Center 
ATTN:  CodeG33. 

T.  Ekrtan 
J.  Copley 

Code  G30.  Guns  and  Munitions 
Division 

Code  G301.  D.  Wilson 
Code  G32,  Gun  Systems  Division 
Code  E23,  Techni^  Library 
Dahlgien.VA  22448-5000 

Commander 

Naval  Underwater  Systems  Center 
Energy  Conversion  Dqpt. 

ATTN:  Technical  Library 
Newport.  RI  02840 

Commander 

Naval  Surface  Warfare  Center 
Indian  Head  Division 
ATTN:  610.  C.  Smith 
6110J.  K.  Rice 
6nOC.  S.  Peters 
Indian  Head.  MD  20640-5035 

Commander 

Naval  WetgxNis  Center 

ATTN;  Code  388.  C.  F.  Price 

Info  Science  Div 

China  Lake.  CA  93555-6001 


1  OSD/SDIOAST 

ATTN:  Dr.  Len  Cavcny 
Pentagon 

Washington.  DC  20301-7100 

1  OLAC  PUTSTL 
ATTN;  D.  ShipleU 
Edwards  AFB,  CA  93523-5000 

1  Commandant 

U.S.  Army  Field  Artillery  School 
ATTN:  STSF-TSM-CN 
Fort  Sill.  OK  73503-5600 

10  Central  Intelligence  Agency 

Office  of  Central  Reference  Dissemination  Branch 
Room  GE47  HQS 
Washington.  DC  20502 

1  Central  Intelligence  Agency 
ATTN:  Joseph  E.  Backofen 
NHB,  Room  5N0I 
Washington.  DC  20505 

5  Director 

Sandia  National  Laboratories 
ATTN:  T.  Hitchcock 

R.  Woodfln 
D.  Benson 

S.  Kempka 
R.  Beasley 

Advanced  Projects  Div  14 
Organization  9123 
Albuquerque,  NM  87185 

2  Director 

Los  Alamos  National  Laboratory 
ATTN:  B.  Kaswhia 
H.  Davis 

Los  Alamos,  NM  87545 
I  Director 

Lawrence  Livermore  National  Laboratory 
ATTN:  M.  S.  L-355,  A.  Buckingham 
P.O.  Box  808 
Livermore,  CA  94550 
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2  Director 

Sandia  National  Laboratories 
Combustion  Research  I^ility 
ATTN:  R.  Armstrong 
S.  Vosen 
Division  8351 

Uvermore,  CA  94SS1-0469 

1  University  of  Dlinr^ 

Dept  of  MecfaTIndust  Engr. 

ATIN:  Professor  Herman  Kiier,  144  MEB 
1206  N.  Green  Sl 
Ufb8na.IL  61801 

1  The  Johns  Hofddns  UniversityA71A 
ATTN:  T.  Christian 
10630  Little  Patuxent  Parkway,  Suite  202 
C<riumbia.MD  21044-3200 

1  Pennsylvania  State  University 
Dept  of  l^dranical  Engr. 

ATTN:  Dr.  K.  Kuo 

312  htohanical  Engineering  Bldg. 

University  Park.  PA  16802 

1  North  Carolina  State  University 

ATTN:  John  G.  GOligan 
60x7909 

1110  Burlington  Engineering  Labs 
Raleigh.  NC  27695-7909 

1  SRI  International 

Propulsion  Sdences  Division 
ATTN:  Technical  Library 
333  Ravenswood  Ave. 

Menlo  Park.  CA  94025 

1  SPARTA 

ATTN:  Dr.  Michael  Holland 
9455  Towne  Center  Dr. 

San  Diego.  CA  92121-1964 

5  FMC  Naval  Systems  Division 

ATTN:  Mr.  G.  Johnson 
Mr.  M.  Seale 
Dr.  A.  Giovanetti 
Mr.  J.  Dyvik 
Dr.  D.  Cook 
4800  East  River  Rd. 

Minneapolis,  MN  55421-1498 
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3  GT  Devices 

ATTN:  Dr,  S.  Goldstein 
Dr.  R.  J.  Greig 
Dr.  N.  Winsor 

5705A  General  Washington  Dr. 

Alexandria.  VA  22312 

3  General  Dynamics  Land  Systems 
ATTN:  Dr.  B.  VanDeusen 
Mr.  F.  Lunsford 
Dr.  M.  Weidner 
P.O.  Box  2074 
Warren.  MI  48090-2074 

2  Science  i^licadons  International  Corporation 
ATTN:  Mr.  N.  Sinha 
Dr.  S.  Dash 

501  Office  Center  Drive.  Suite  420 
Fbit  Wadiington,  PA  19034-3211 

2  AUiam  Techsystems,  Inc. 

ATTN:  R.E.  Tompkins 
J.  Kemiedy 
MN38-3300 

10400  Yellow  Circle  Dr. 

Minnetonka.  MN  55343 

2  Olin  Ordnance 

ATTN:  V.  McDonald,  Library 
Hugh  McEiroy 
P.O.  Box  222 
Sl  Marks,  FL  32355 

1  Gough  Associates.  Inc. 

ATTN:  P.  S.  Gough 
1048  South  Sl 

Portsmouth,  NH  03801-5423 

1  Physics  International  Library 
ATTN:  H.  Wayne  Wamider 
P.O.  Box  5010 

San  Levidro.  CA  94577-0599 

2  Rockwell  International 
Rocketdyne  Division 
ATTN:  BA08. 

J.  E.  Flarugan 
J.  Gray 
6633  Canoga  Ave. 

Canoga  Park,  CA  91304 
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Aberdeen  Proving  Ground 

4  Cdr.  USACSTA 
ATTN;  S.  Walton 
G.  Rice 
D.  Lacey 
C.  Herud 


2  Science  Applications  International  Corporation 
ATTN:  J.  Batteh 

L.  Thornhill 
1S19  Johnson  Ferry  Rd 
Suite  300 

Marietta.  GA  30062-6438 

1  Eli  Freedman  &  Associate 

ATTN:  E.  Freedman 
2411  Diana  Rd. 

Baltimtue.  MD  21209 

1  Veriiay  Technology.  Inc. 

4845  Millersp^Kt  Hwy. 

P.O.  Box  305 

Faa  Amherst.  NY  14051*0305 

1  Battelle 

ATTN;  TACIEC  Library.  J.  N.  Huggins 
505  King  Ave. 

Columbus.  OH  43201-2693 

2  California  Institute  of  Technology 
Jet  Propulsion  Laboratory 
ATTN:  L.  D.  Strand.  MS  125-224 

D.eUot 

4800  Oak  Grove  Dr. 

Pasadena.  CA  91109 

1  General  Electric  Co. 

Defense  Systems  Division 
ATTN:  Dr.  J.  Mandzy 
MaU  Drop  43-220 
100  Plastics  Ave. 

Pittsfield.  MA  01201 

1  Dept  Electrical  Engineering 

ATTN:  Dr.  W.  J.  Sargeant, 

James  Clark  Maxwell  Professor 
Bonner  Hall  -  Room  312 
Buffalo.  NY  14260 
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2  Princeton  Combustion  Research 
Laboratories 

ATTN:  M.  Summerfield 
N.  Messina 

Princeton  Corporate  Plaza 
11  Deerpark  Drive 
Bldg.  IV.  Suite  119 
Monmouth  Junction.  NJ  08852 
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2  RARDE 

GS2  Division 
Building  R31 
ATTN:  Dr.  C.  Woodley 
Dr.  G.  Cook 
Fort  Halstead 

Sevenoaks,  Kent  TN14  7BP 
En^and 


iNTEMnONALLY  LEFT  BLANK. 
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USER  evaluation  SHEET/CHANGE  OF  ADDRESS 

This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes.  Your 
comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  ARL  Report  Number  ARL-TR-6 3 _ Date  of  Report  February  1993 _ 


2.  Date  Report  Received  _ _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for 

which  the  report  will  be  used.)  _ _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  dau,  procedure,  source  of 
ideas,  etc.)  _ _  - _ 


S.  Has  the  information  in  this  repon  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved, 
operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  General  Conunents.  What  do  you  think  should  be  changed  to  improve  future  rqx>rts?  (Indicate 
changes  to  organization,  technical  content,  format,  etc.) _ 


Organization 

CURRENT  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 

City,  Sute,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address 
above  and  the  Old  or  Incorrect  address  below. 


Organization 

OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 

City.  State.  Zip  Code 

(Remove  this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  mail.) 


DEPAfTTMCNT  OF  THE  ARMY 


OFROAL  BUSINESS 


BUSINESS  REPLY  IVL^IL 

BRST  CUBS  mm  Mb  DOm,  tfS.  MD 

Postage  «ill  bt  paid  by  addiessec 


Director 

U.S.  Army  Research  Laboratory 
ATTN:  AMSRL-OP-CI-B  (Tech  Ub) 
Aberdeen  Proving  Ground,  MD  21005-5066 


